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ABSTRACT 


Seafloor  morphology  in  the  surf  zone  of  a  sandy  ocean 
beach  was  measured  nearly  continuously  for  45  days  with  a 
1.4  X  1.4  m  coherent  array  of  7  sonar  altimeters.  Migrating 
megaripples,  seafloor  bedforms  with  amplitudes  of  0(10-50 
cm)  and  lengths  of  0(1-5  m)  ,  were  observed  in  about  2  m 
water  depth  in  the  trough  between  a  sand  bar  and  the 
shoreline  for  a  wide  range  of  wave  and  current  conditions . 
Megaripple  migration  speed  and  direction  are  estimated  from 
the  array  data  using  cross-correlations  between  seafloor 
elevation  time  series  observed  along  the  cross-  and 
alongshore  array  legs.  Megaripples  were  shown  to  be  aligned 
in  a  direction  that  maximized  the  gross  sediment  transport 
normal  to  the  bedform  crest  (Rubin  and  Hunter,  1987; 
Gallagher,  et  al .  1998).  It  is  hypothesized  that  megaripple 
migration  rate  is  related  to  the  net  transport  in  the 
direction  of  bedform  alignment.  The  speed  of  megaripple 
migration  is  compared  with  the  magnitude  of  the  velocity 
field  normal  to  the  bedform  crest  in  the  direction  of  the 
mean,  wave  orbital,  and  resultant  velocities. 
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I .  INTRODUCTION 


Megaripples  are  bedforms  with  heights  of  10  -  50  cm  and 
lengths  of  1  -  5  m.  These  features  are  important  because 
they  affect  bedload  transport,  sediment  suspension,  boundary 
layer  turbulence  generation,  energy  dissipation  of  waves, 
and  currents  via  form  drag.  Few  observations  of  bedforms  in 
the  nearshore  exist  because  the  surf  zone  is  a  difficult 
environment  in  which  to  make  measurements.  Clifton  et  al. 
(1971)  and  Clifton  (1976)  made  measurements  of  bedforms  on 
beaches  in  Oregon  using  SCUBA  during  mild  conditions.  They 
observed  megaripples  (Figure  1)  both  inside  the  surf  zone 
(i.e.  in  the  bar-trough)  and  outside  the  breaking  waves 
(Figure  2).  Clifton  et  al .  (1971)  determined  lunate 
megaripples  migrated  towards  the  shore  at  an  average  rate  of 
30  cm/hr  during  low-energy  conditions.  Clifton  (1976) 
proposed  the  onshore  migration  was  a  result  of  sediment 

transport  associated  with  asymmetries  in  wave  velocities. 

\ 

They  concluded  that  these  structures  are  common  in  medium- 
to  coarse-grained  sand  under  conditions  of  intense 
asymmetric  flow  generated  by  long  period  waves  and  are 
common  in  2  -  4  m  of  water.  Hay  and  Wilson  (1994)  observed 
similar  crescent  shaped  bedforms  using  a  fixed  rotating 
sidescan  sonar.  Hay  and  Wilson  (1994)  and  Hay  and  Bowen  (in 
press)  observed  migration  rates  between  0.5  to  3  cm/hr 
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during  both  low  and  high  wave  energy  conditions .  During  a 
storm  Hay  and  Bowen  (in  press)  found  a  strong  relationship 
between  alongshore  current  speed  (20  to  80  cm/s)  and 
alongshore  megaripple  migration  rate,  0(3  cm/hr) .  Gallagher 
et  al.  (1998b)  compared  megaripple  migration  direction  with 
the  direction  of  waves,  mean  currents,  and  their  vector  sum 
as  well  as  with  a  model  for  desert  dune  orientation  (Rubin 
and  Hunter,  1987)  .  They  found  that  the  model,  which  is 
based  on  instantaneous  flow  vectors,  was  a  better  predictor 
of  migration  direction  than  any  time  averaged  flow  parameter 
(mean,  root  mean  square  (rms),  and  their  vector  sum) . 

If  migration  of  bedforms  is  an  important  mechanism  for 
bedload  sediment  transport,  then  the  understanding  of 
megaripple  migration  in  oscillatory  flow  is  fundamental  to 
understanding  nearshore  hydrodynamic  and  sediment  processes. 
When  they  migrate,  they  may  be  an  important  mechanism  for 
sediment  transport.  The  relationship  between  the  migration 
rate  of  megaripples  and  wave  and  current  forcing  is  studied 
in  this  thesis. 
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II .  DUCK94  EXPERIMENT 


The  measurements  were  obtained  during  DUCK94,  between 
01  September  1994  through  15  October  1994  at  the  U.S.  Army- 
Corps  of  Engineers  Field  Research  Facility  (FRF)  in  Duck, 
North  Carolina.  The  DUCK94  field  experiment  was  a  pilot 
study,  which  represented  a  first  attempt  to  make 
quantitative'  multi-scale  measurements  of  nearshore 
morphology  (Thornton  et  al.,  1998).  The  FRF  is  located  on 
the  Outer  Banks,  a  barrier  island  formation  with  no  major 
coastal  structures  that  can  interfere  with  nearshore  flows. 
The  mean  amplitude  of  the  semidiurnal  tide  is  approximately 
0.5  m.  Fine  quartz  sediments  within  the  surf  zone  are 
fairly  well  sorted  with  a  mean  grain  size  of  0.2  mm  while 
sediments  on  the  foreshore  are  poorly  sorted  with  larger 
mean  grain  sizes  usually  greater  than  0.4  mm. 

Observations  of  migrating  megaripples  were  made  in  the 
surf  zone  using  an  array  of  seven  1  MHz  downward- looking 
automatic  gain  control  (AGC)  sonar  altimeters  (Figure  3)  . 
Sonar  altimeters  estimate  the  location  of  the  seafloor  by 
measuring  the  time  for  an  acoustic  pulse  to  travel  to  the 
seafloor  and  back  and  then  converting  the  time  to  a  distance 
using  the  theoretical  speed  of  sound.  The  AGC  can 
accurately  estimate  the  seafloor  location  to  within  ±  3  cm 
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even  in  the  surf  zone  where  energetic  breakers  and  strong 
currents  produced  bubbles  and  dense  sediment  suspension 
(Gallagher  et  al.,  1996). 

The  array  was  located  approximately  50  m  from  the 
shoreline  in  1.5  to  2  m  of  water  depth  (Figure  4)  in  the 
trough  between  the  sand  bar  and  the  beach  (Figure  3a) .  Wave 
height  ranged  from  0.5  to  4  m  with  three  large  storms  (02  - 
05  September,  02  -  04  October,  and  10  -  20  October) 

occurring  during  the  experiment  (Figures  5  and  6)  .  An 
electromagnetic  current  meter,  which  measured  cross-  and 
alongshore  components  of  the  fluid  velocity  at  2  Hz  was  co¬ 
located  with  the  array  and  was  about  0.5  m  above  the  bed 
(Figure  3b)  .  During  storms,  both  wave  orbital  velocities 
(indicated  by  Urms)  and  steady  currents  {Umeon,  Vmean)  were 
large,  0(1  m/s)  (Figure  7).  During  periods  of  low  energy 
waves  (wave  height  <  1  m)  ,  currents  were  negligible  and  RMS 
currents  were  small,  0(25  cm/s) .  Megaripples  were  observed 
to  exist  about  60%  of  the  time  and  under  a  variety  of 
conditions . 
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III.  THEORY  AND  METHODOLOGY 


Traditionally,  sediment  bedforms  in  directionally 
varying  flow  have  been  classified  as  transverse, 
longitudinal,  or  oblique  by  their  alignment  with  the 
resultant  transport  vector  (normal,  parallel,  or  neither) . 
Rubin  and  Hunter  (1987)  proposed  that  gross  transport  should 
be  considered  since  all  transport  can  build  bedforms  even  if 
the  net  or  resultant  is  zero  (i.e.  purely  oscillatory  flow) . 
They  hypothesized  that  bedforms  are  aligned  such  that  the 
gross  sediment  transport  normal  to  their  crests  is 
maximized. 

To  predict  the  orientation  of  the  bedforms,  the 
absolute  value  of  the  components  of  the  two  transporting 
flows  normal  to  the  bedform  crests  are  added, 

T  =  D\sin a\  +  S\sin(y  -cc]  ( 1 ) 

where  D  and  S  are  the  dominant  and  subordinate  vectors  of 
the  two  flows,  a  is  the  angle  between  D  and  the  bedform 
crest,  and  y  is  the  divergence  angle  between  D  and  S  (Figure 
8) .  The  first  derivative  is  applied  to  find  the  value  of  a 
which  maximizes  T, 
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tana  = 


± 


D\  , 

— r  +  \cos 


y\ 


sin  y| 


(2) 


Rubin  and  Hunter  (1987)  tested  the  alignment  of  subaerial 
bedforms  in  bi-directional  flow  by  using  a  rotating  bed  in  a 
unidirectional  flow.  They  found  good  agreement  between 
their  hypothesis  and  the  rotating  bed  experiment.  The 
results  showed  that  for  a  given  pair  of  flow  vectors,  D  and 
S,  the  bedforms  aligned  such  that  the  gross  bedform  normal 

transport  was  maximized.  Rubin  and  Ikeda  (1990)  studied  the 
alignment  of  subaqueous  dunes  in  a  laboratory  experiment  and 
similarly  found,  "The  bedforms  had  the  orientation  that  was 
subject  to  more  transport  crossing  the  bedform  Crestline 
than  any  other  possible  bedform  orientation."  They  also 
learned  two  important  requirements  of  this  model:  1)  the 
transporting  vectors  must  be  temporally  distinct,  and  2) 
they  must  be  relatively  short  enough  in  duration  that  the 
bedforms  do  not  come  into  equilibrium  with  either 
transporting  flow. 

Gallagher  et  al .  (1998b)  applied  this  theory  to  field 
observations  made  in  the  surf  zone  during  the  DUCK94 
experiment.  Observed  megaripple  migration  rate  was 
calculated  quantitatively  by  estimating  the  time-lag  between 
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altimeter  pairs  from  the  maximum  cross-correlation.  A 
sample  of  the  time  series  from  three  sonars  in  the  cross¬ 
shore  leg  of  the  array  are  shown  in  Figure  9.  The  three 
panels  are  positioned  relative  to  the  actual  spacing  (where 
the  top  panel  is  from  the  most  onshore  altimeter)  and  the 
line  drawn  between  similar  features  illustrates  a  30  cm/hr 
migration  rate.  Overlapping  48-hour-long  time  series  (e.g. 
Figure  9)  were  used  to  estimate  megaripple  migration  speeds 
and  directions  for  the  two  month  long  experiment  (when 
megaripples  existed) . 

To  apply  Rubin  and  Hunter  (1987)  to  nearshore  forcing, 
D  and  S  need  to  be  properly  represented.  In  the  nearshore, 
mean  and  oscillatory  currents  are  traditionally  separated 
and  treated  as  different  forcing  mechanisms  for  sediment 
transport.  However,  they  are  not  temporally  distinct,  and 
thus  can  not  be  used  to  represent  D  and  S.  Gallagher  et 

al .  (1998b)  used  the  onshore  wave  surge  (plus  any  coincident 
mean  flow)  and  the  offshore  wave  surge  (plus  any  coincident 
mean  flow)  to  represent  D  and  S  and  to  predict  orientation 

of  bedf orms .  This  representation  gives  temporally  distinct 
vectors.  In  addition,  the  time  required  for  megaripples  to 

ibrate ,  O(hrs)  ,  is  long  compared  with  wave  velocities, 
0 (secs) ,  thus  the  model  is  valid  for  use  in  the  surf  zone. 

To  parameterize  the  onshore  and  offshore  surge  of  the 
waves  from  measured  currents  sampled  at  2  Hz,  instantaneous 
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transport  was  approximated  by  the  magnitude  of  the 
instantaneous  velocity  cubed 

Q  ^(f/2+y2)  (3) 

(after  common  transport  models,  e.g.  Bailard,  1981) .  These 
instantaneous  transport  vectors  were  sorted  into  5°  bins 
then  siammed  over  a  3 -hour  period  given  a  directional 
distribution  of  cumulative  transport  (Figure  10)  .  Using 
these  cumulative  transport  vectors,  bedform  orientation  was 
found  by  numerically  solving  for  bedform  orientation,  a, 
for  which  the  gross  bedform  normal  (GBN)  transport  was 
maximized 


72 

Tgbn  =  X  Q-^sinifc  —  7,)| 

;=i 


(4) 


where  i  represents  each  5°  bin.  The  3 -hour  estimates  were 

then  averaged  to  give  overlapping  48-hour  estimates  of 
bedform  orientation  for  comparison  with  observed  megaripple 
migration  direction. 

The  observed  megaripple  migration  direction,  Qunp,  are 
compared  with  the  predicted  direction,  6rh,  in  Figure  11 
where 
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6rh  -  90°  a  . 


(5) 


Gallagher  et  al .  (1998b)  found  that  the  oriention  of 

nearshore  megaripples  is  explained  by  Rubin  and  Hunter's 
model  of  subaerial  dune  migration  direction. 

Here,  it  is  hypothesized  that  megaripple  migration  rate 
is  related  to  the  net  transport  in  the  direction  of  bedform 
alignment.  This  net  bedform  normal  transport  is  calculated 
by  summing  the  directional  bin  components,  which  are  normal 
to  the  predicted  Grh  such  that 


72  ,  . 

Net  T gbn  =  X  Qi  sin\a  -  yj  ( 6 ) 

i=i 

/ 

where  i  represents  each  5°  bin.  Note  that  it  is  now 

possible  to  obtain  positive  and  negative  migration  rates 
since  we  include  the  directionality  of  the  normal  component. 
Again,  3 -hour  estimates  of  megaripple  migration  rate  were 
averaged  to  give  overlapping  48-hour  estimates  for 
comparison  with  observed  megaripple  migration  rates . 
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IV.  ANALYSIS 


Measured  velocities  are  used  to  parameterize  the  flow 
field  for  comparison  with  observed  megaripple  migration 
rates.  The  measured  cross-  (U)  and  alongshore  ( V )  fluid 

velocities  were  separated  into  mean  and  wave  (oscillatory) 
1-hour  averages 


(7) 

(8) 

where  the  bar  and  the  tilde  represent  the  mean  and  wave 
velocities  respectively.  Next,  hourly  skewness-weighted  rms 
wave  velocities  were  calculated 

(9) 

(10) 


y  ~  _ _ 

skew  /~?\ 
(v  ) 


U  -  u  +  u 
V  =  v  +  v 


where  (  )  denotes  time  averaging.  Hourly  mean  and  skewness- 

weighted  rms  wave  velocities  were  then  averaged  for  three 
hours 
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(11) 


«3»V3  =  (m),(v) 

=  (uskew),(vskew)  .  (12) 


Magnitudes  for  3 -hour  mean  velocities,  skewness- 
weighted  rms  wave  velocities,  and  the  vector  sum  of  the  mean 
and  wave  velocities  were  calculated 


TJ  3  —  +Vj 
U3  =  V«32  +  V32 


C/3  =  «/(h3  +  m3)2  +  (v3  +  v3  . 


(13) 

(14) 

(15) 


Finally,  48-hour  magnitudes  of  the  mean,  the  skewness- 
weighted  rms  wave,  and  the  vector  sum  were  calculated  for  16 
sequential  3 -hour  magnitudes 


vtt  =  (v,) 

(16) 

ro 

ii 

00 

to 

(17) 

<5 

-p^ 

oo 

II 

"cf 

U) 

(18) 

in  order  to  compare  with  the  48-hour  observed  megaripple 
migration  rates. 
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V.  RESULTS  AND  CONCLUSIONS 


Observed  megaripple  migration  speeds  are  plotted  in 

Figure  12  verses  the  cube  of  the  mean  ( f/48 )  ,  skewness- 

weighted  rms  wave  (U4S)  ,  and  vector  sum  (t/48)  velocities  in 
addition  to  the  predicted  net  gross  bedform  normal 

transport,  Net  Tgbn ■  Results  were  not  sensitive  to  the 
exponent  (2,  3,  4,  or  6)  in  the  model  used  to  represent 

instantaneous  sediment  transport  ( Q  in  Equation  3  and  Qi  in 
Figure  10)  .  The  correlation  coefficient  values,  r,  are 

similar  for  U4i,  £/48 ,  and  Net  Tgbn >  and  all  the  r  values  are 

statistically  significant  at  the  95%  confidence  level. 

Therefore,  results  do  not  suggest  that  one  transport 
formulation  is  best. 

There  are  roughly  two  populations  of  data  points.  One 
corresponds  to  a  current-dominated  period  where  megaripple 
speeds  are  greater  than  about  60  cm/hr,  while  the  other 
period  is  wave-dominated  with  megaripple  speeds  less  than 
about  60  cm/hr.  In  panels  (a)  ,  (c)  ,  and  (d)  ,  the  current- 

dominated  points  all  correspond  to  large  transport 
estimates,  so  megaripple  speeds  could  be  dependent,  at  least 
crudely,  on  any  of  these  three  transport  parameters .  In 
general,  these  observations  suggest  that  high  megaripple 
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speeds  correspond  to  high  current  velocities  when  sediment 
transport  is  large  .and  low  speeds  correspond  to  low  current 
velocities  when  less  sand  is  in  motion.  This  result  is  in 
general  agreement  with  the  results  of  Hay  and  Bowen  (in 
press)  who  also  observed  migrating  megaripples  during 
DUCK94.  They  focused  on  a  three  day-long,  high-energy  event 
during  which  they  found  that  megaripple  migration  rates  were 
a  strong  function  of  the  observed  steady  currents. 

For  low-energy  periods  when  oscillatory  wave  velocities 
are  dominant  and  megaripple  migration  is  slow,  there  is  no 
strong  relationship  between  megaripple  speed  and  t/48  (Figure 

12b) .  Improved  models  for  sediment  transport,  which  perhaps 
include  acceleration  (Clifton,  1976;  Gallagher  et  al . , 
1998a),  might  improve  this  relationship. 

These  results  suggest  that  during  periods  of  high- 
energy  waves,  megaripple  migration  rates  are  higher  than 
during  periods  of  low-energy  waves.  Thus,  the  hypothesis 
that  megaripple  migration  rate  is  related  to  transport  is 
supported  by  the  analysis  of  this  data  set.  The  hypothesis 
that  net  transport  in  the  direction  of  bedform  alignment 
following  Rubin  and  Hunter  (1987)  and  Gallagher  et  al . 
(1998b)  would  be  the  best  predictor  of  megaripple  speed  is 
not  supported.  However,  the  prediction  of  sediment 
transport  via  megaripple  migration  requires  a  vector,  that 
is,  both  direction  and  magnitude.  Thus,  the  model  of  Rubin 
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and  Hunter  (1987)  is  a  better  tool  for  predicting  bedform 
migration.  More  study  in  this  area  is  required  in  order  to 
better  understand  bedform  processes  and  sediment  transport 
in  the  nearshore. 
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APPENDIX:  FIGURES 


17 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


18 


Vertical  Exaggeration  2X 


Figure  1.  Plan  view  and  profile  across  lunate  megarippl 
(Clifton  et  al . .  1911). 
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Figure  2.  Megaripples  have  been  observed  in  the  outer  rough 
facies  and  the  inner  rough  facies  of  a  high-energy  nearshore 
(Clifton  et  al . ,  1911). 
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Figure  3.  Location  of  the  altimeter  array.  (a)  X  shows  the 
location  of  the  array  in  the  trough  between  the  sandbar  and 
the  beach.  The  array  was  approximately  50  m  from  ‘  the 
shoreline.  Profiles  were  from  01  SEP  1994  (solid  curve)  and 
15  OCT  1994  (dashed  curve)  .  (b)  The  plan  view  of  the 

altimeter  array  shows  the  location  of  the  AGO  altimeters 
(circles),  the  electromagnetic  current  meter  (square),  and 
the  three  altimeters  used  in  Figure  (9)  time  series 
(asterisks). 
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Figure  5 .  Climatology  during  September  for  the  DUCK94 
experiment.  H ^  is  the  spectral  estimate  of  the  significant 
wave  height.  Tp  is  the  period  of  the  peak  frequency 
measured  in  8  m  depth. 
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Figure  6 .  Climatology  during  October  for  the  DUCK94 
experiment.  Hmo  is  the  spectral  estimate  of  the  significant 
wave  height.  Tp  is  the  period  of  the  peak  frequency 
measured  in  8  m  depth. 
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Figure  7 .  Three  hour  averages  of  the  instantaneous  fluid 
velocity  components  from  01  SEP  1994  to  15  OCT  1994. 
Positive  values  indicate  onshore  or  southerly  flow  and 
negative  values  indicate  offshore  or  northerly  flow. 
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Figure  8.  Transport  vector  example  (Gallagher  et  al . , 
1998b)  where  7  is  the  angle  between  the  dominate,  D,  and 
subordinate,  S,  transport  vectors.  cc  is  the  angle  between 
D  and  the  bedform  crest  (double  line) .  R  is  the  vector  sum 


of  D  and  S. 
transport . 


is  the  gross  bedform  normal 


Figure  9.  Migrating  megaripples  are  seen  in  this  time 
series  example  from  three  altimeters  (asterisks  in  Figure 
3b)  where  (a)  is  from  the  most  onshore  altimeter.  The 
slopes  of  the  lines  connecting  the  troughs  of  the  bedforms 
illustrate  onshore  migration  on  0(30  cm/hr)  (Gallagher  et 
al.,  1998b). 


Figure  10.  Example  of  predicted  bedform  orientation  where 
the  magnitude  and  direction  of  each  cumulative  transport 

vector  are  given  by  Qi  and  y.,  respectively,  a  is  the 
orientation  of  the  bedform  crest  (double  line)  such  that  the 
gross  bedform  normal  transport,  T,  is  maximized.  6 m  =  90°  - 

a.  The  direction  of  the  resultant  vector,  R,  is  included 
for  comparison  (Gallagher  et  al . ,  1998b). 
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Figure  11.  Predicted  megaripple  migration  direction,  Q^, 
versus  observed  megaripple  migration  direction,  Q  where 

0°  is  onshore  and  90°  is  to  the  south  and  r  is  the 
correlation  coefficient  (Gallagher  et  al . ,  1998b). 
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